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Abstract. 

In this paper, we discuss why and how the Planck mission, originally designed 
and proposed for mapping intensity fluctuations, has been revised for polarisation 
measurement capability as well. 
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1. Introduction 

After the great success of Cosmic Microwave Background (CMB) anis- 
otropics as a data goldmine for constraining cosmological models, a 
large fraction of the CMB experimental community now turns towards 
building experiments for measuring the small fraction of polarised emis- 
sion in the CMB. Motivations for doing so are many, but the experi- 
mental effort quite challenging. 

Polarisation signals are much weaker than temperature anisotropics. 
The amplitude of polarisation fluctuations from scalar modes is about 
one tenth of that of CMB anisotropics on small scales, and much lower 
(relatively) on large scales. The most interesting polarisation signature, 
that of gravity waves from inflation, is even smaller by yet an order of 
magnitude, possibly more. Hence, very sensitive instruments are needed 
to measure the polarisation signals of the CMB. 

The Planck mission, to be launched by ESA in 2007, is the third gen- 
eration CMB mission after COBE and WMAP. In this paper, we review 
the motivations for CMB polarisation observations and the design and 
expected outcome of the Planck mission for polarisation measurements. 



2. CMB polarisation: Theoretical review 

The standard paradigm for the generation of CMB anisotropics is now 
well understood. Initial perturbations, generated during an inflationary 
epoch in the very early Universe, evolve under the effect of gravity, 
which causes over-densities to collapse, opposed by radiation pres- 
sure, which prevents structures from collapsing as long as radiation 
(mainly photons) is tightly coupled to matter. These opposing forces 
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generate acoustic oscillations in the photon-baryon fluid prior to de- 
coupling, giving rise to peaks and troughs in the anisotropy spatial 
power spectrum. 

2.1. The generation of CMB polarisation 

Temperature anisotropics on the sky are essentially generated through 
the combination of three effects. First, regions of the Universe in- 
trinsically hotter at last scattering emit more and hotter photons. 
Second, photons coming out of overdensities are redshifted by gravi- 
tational effect, while those coming from underdensities are blueshifted. 
Finally, fluid motions at last scattering generate an additional redshift 
or blueshift by Doppler effect. 

On large scales, causal processes have had no time to change much 
the primordial fluctuations. Hence, the CMB anisotropics reflect essen- 
tially the original perturbations (with some modifications due solely to 
the Integrated Sachs- Wolfe effect on large scales). On sizes smaller than 
the horizon, however, original perturbations have been modified by a 
transfert function which depends on the physics governing the acoustic 
oscillations. Hence, the relative height of the acoustic peaks depend 
on the mass/energy density content of the Universe (of baryons which 
are the main source of inertia, of photons which are the main source 
of radiation pressure, of dark matter which contributes most of the 
depth of the potential wells during the matter domination epoch) and 
on the strength of the damping term due to the expansion (and hence 
on the Hubble constant). Finally, the translation from physical sizes 
to angular sizes depends on the geometry, and hence on cosmological 
parameters as the curvature and the cosmological constant. 

To some extent, the generation of polarisation is simpler than that 
of temperature anisotropics. Polarisation can be generated solely at last 
scattering by the dependence of the Thomson cross section upon the the 
polarisation states of the incoming and outgoing photons, aT oc le^.e^p. 
Hence, a quadrupole in the radiation intensity at some point at last 
scattering generates polarisation in a direction perpendicular to that 
of the maximum of the quadrupole tangentially to the last scattering 
surface (Hu & White, 1997). 

Contrarily to CMB temperature fluctuations which arise from a su- 
perposition of monopoles, dipoles (from Doppler effect) and quadrupoles 
at last scattering, CMB polarisation arises from quadrupoles alone. 
Hence, polarisation is an essential tool to separate between sources of 
anisotropics (monopoles, dipoles, quadrupoles), and lift some of the 
possible degeneracies among cosmological parameter sets as compared 
with constraints obtained from CMB anisotropics alone. 
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2.2. Polarisation description 



2.2.1. Stokes parameters 

The polarisation state of an incoming transverse electromagnetic wave 
is fully described in a given reference frame by four parameters /, Q, 
[/, and V defined as follows: 



^.\') + {\Ey\') 



Q = {\E,\^) - {\Ey\^) 

U = {E^eD + {EyEl) 
V = i{{E,El) - {EyEl)) 



(1) 
(2) 
(3) 
(4) 



/, Q, U ^ and V are known as the Stokes parameters of the radiation. 
/ gives the total intensity, Q and U describe the linearly polarised part, 
and V the circularly polarised part of the radiation. 

The measured Stokes parameters depend on the reference system 
(X, Y). In a frame (X', Y') obtained from (X, Y) by a rotation by an 
angle iIj^ the Stokes parameters Q' ^ U' ^ and V are obtained from /, 
Q, f7, and V by the following transformation: 
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Hence, Q and U are defined in a frame-dependent way. At each point 
of the sphere, typically, directions along spherical coordinates 9 and 
(j) can be used as reference - although other choices can be made. As 
a consequence, changing frames from ecliptic to galactic coordinates 
(and redefining correspondingly coordinates 6 and (j) and polarisation 
reference axes at each point) requires computing in each pixel of the 
new map a position dependent linear combination of Q and U in the 
corresponding pixel of the original map. More generally, this will be the 
case in any map transform which does not preserve reference axes direc- 
tions. Local projections onto small square maps for instance, especially 
near the poles where there is a coordinate singularity, involve the same 
kind of difficulties. Hence, the reprojection and/or repixelisation is even 
more technical for polarisation maps than for temperature maps. 

2.2.2. Spherical harmonics expansion 

It has now become traditional to expand CMB temperature fluctuations 
on a basis of spherical harmonics on the sky. Similarly, it is desirable 
to expand polarisation fields. A method for doing so has been proposed 
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by Seljak and Zaldarriaga (1997), and independently by Kamionkowski, 
Kosowski and Stebbins (1997). 

The fact that Q and U depend on the coordinate system at each 
point of the sphere is a comphcation for a meaningful expansion. How- 
ever, combinations Q±iU (known as the Kuscer parameters (Kuscer & 
Ribaric, 1959)) transform, under a rotation of coordinates by an angle 
Q^, as Q zb iU — > e^'^^^{Q ± iU). Hence, Q ±iU are eigenmodes of the 
"rotation by angle a" operator on the plane, associated to eigenvalues 
of e=F22Q; They are spin ±2 quantities, and can be decomposed on the 
sphere on the basis of spin 2 spherical harmonics ±2Y£m as: 

Q{9, 0) ± iU{9, 0) = ^ a^2/m ±2Yim{0, 0) (5) 

i,m 

The E and B polarisation fields can be defined as 

E = y2ag,Yem{e,(f>) (6) 



and 

where 
and 



B = Taf^YUO,'l>) (7) 



(^£m = ■ 7^ — (8) 



«£m = 2 

The E and B fields are respectively a scalar and a pseudo-scalar field 
on the sphere (i.e. fields with respectively even and odd parity). Note 
that the definition of af^ and af^ given in equations 8 and 9 is just 
a particular choice for obtaining a scalar and pseudo-scalar field from 
a2/rn and a_2/m- Any other choice af^ — > a{£)af^ would be just as 
good, but we stick here to the definition of Seljak and Zaldarriaga - 
which is also the one used in the CMBfast software. 

The intensity and polarisation of the sky is fully described by the 
set I^E^By . In turn, the spatial power spectrum of the intensity (tra- 
ditionally expressed as a temperature fluctuation in units of /xK, and 
noted T) and polarisation field (in the same units) of the CMB is given 
by: 



( Cj^ Cj^ o\ 
Oj^ Cf^ 
Cf ^ 
0/ 
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In this spatial power spectrum, all V terms vanish because no V- 

type polarisation is expected in the CMB, and all cross terms involving 

B vanish for parity reasons. Hence, the intensity and polarisation power 

spectrum of the CMB is fully described, for each by four quantities, 
riTT r'TE r'EE c'Bb 

For a given cosmological model, it is possible to predict all of these 
terms of the matrix Cf^, which we refer to as the multivariate spatial 
power spectrum of CMB temperature and polarisation. 

2.3. Science with CMB polarisation 

It is clear that polarisation and temperature anisotropics together carry 
more information than temperature anisotropics alone. In this section, 
we discuss the main scientific case for CMB polarisation measurements. 

2.3.1. Reionisation 

As polarisation is generated only at last scattering, it probes last scat- 
tering in a more direct way than anisotropics alone. In particular, 
polarisation provides a strong test for the ionisation history, as a sig- 
nificant reionisation optical depth will result in polarisation generation 
closer to the observer than the last scattering surface, at a later time. 
Large scale polarisation patterns are generated in this way, detectable 
as low-^ bumps in the polarisation power spectra. Recently, the WMAP 
team has announced a surprisingly early reionisation at a redshift of 
about 20 as adjusted on the best fit optical depth r :^ 0.17 (Spergel et 
al., 2003). 

If confirmed, this result also has an impact on the trade-off between 
sky coverage and sensitivity per pixel for the detection of polarisation 
B-modes, which may provide a case for spaceborne observations rather 
than ground based. As argued by Lewis, Challinor and Turok (2002), 
in absence of reionisation, a deep observation of the polarisation in a 
small few degree by few degree patch of the sky (feasible from ground) 
would be more optimal to detect inflationary B modes. With early 
reionisation, there is a large bump at very low I which can only be 
measured accurately with large sky coverage, and hence from space. 

2.3.2. Cosmological parameters 

Polarisation power spectra carry also information about cosmological 
parameters which is partly redundant and partly complementary to 
that imprinted on CMB temperature anisotropics alone. 

Redundancy in the information comes from the fact that tempera- 
ture and polarisation anisotropics are generated from the same set of 
initial perturbation evolved until last scattering by the same physical 
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processes. Hence, the pattern of peaks and troughs is connected to 
cosmological parameters about as unambiguously in polarisation power 
spectra as in the temperature power spectrum. Complementarity comes 
from the fact that the dependence of each of the spectra (TT, TE^ EE^ 
BB) on cosmological parameters suffers from different degeneracies. 
The joint information from all spectra hence permits to lift some of 
them. In addition, it provides consistency tests on the model. 

Of course, as long as errors on parameter determination originate 
essentially from instrumental noise, the accuracy of the determination 
of polarisation power spectra is significantly worse than that of the 
temperature power spectrum, and most of the cosmological information 
comes from the TT power spectrum. The other extreme is the situation 
in which errors originate essentially from cosmic variance. Then, all 
spectra carry about as much information content, although not on the 
same parameters, as discussed below. 

2.3.3. Constraining inflation 

One of the most interesting aspect of CMB polarisation is the possibility 
to access directly some of the inflationary parameters. 

There is more and more compelling evidence that initial pertur- 
bations indeed occurred during an early phase of rapid accelerated 
expansion. The observation of several acoustic peaks in the CMB spa- 
tial power spectrum, the Gaussianity of the perturbations, and the 
spatial flatness of the Universe, all of which are strong predictions of 
inflation, are now well established. 

Still, there are at present many inflationary models, none of which 
is fully satisfactory yet (in particular because of the lack of insight on 
the nature of the inflationary field), with too few observational means 
of discriminating among them. It turns out that CMB polarisation 
provides strong constraints on a large class of inflationary scenarios 
- and can possibly permit one to measure some of the fundamental 
parameters of the inflationary field. 

Most inflationary models indeed predict that initial perturbations 
comprise a contribution from tensor modes (gravity waves) in addi- 
tion to scalar modes (density perturbations). Tensor and scalar modes 
contribute to the anisotropy spectrum on large scales. 

For a large class of inflationary models (slow-roll models), inflation 
is described by three main parameters (Liddle & Lyth, 1992) 



^ 167r 



V 

v_ 
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r = T/S 

where V{(j)) is a potential for the inflaton field 0, and T and S are 
the amplitudes of the tensor and scalar perturbations respectively. The 
tensor to scalar ratio, r = T/S^ is expected to be of about 0.1 at most 
and depends on the inflationary model (it is connected to the energy 
scale of the inflation). The scalar and tensor spectral indexes, ns and 
nT, are obtained from e and rj as 

= 1 — 6e + 2?7 

nt = -2e 

Hence, for this class of model, measuring the tensor to scalar ratio 
T/S and the spectral indices ns and constrains the inflationary 
potential V{(f)) and the physics of inflation. 

While ns can be measured from temperature anisotropics, the mea- 
surement of the contribution of tensor modes to the TT power spectrum 
is limited by cosmic variance — essentially because tensor modes con- 
tribute power only on large scales {£ < 100 or so) as they decay upon 
entering the horizon. Small values of T/S^ therefore, are not measurable 
on the TT power spectrum. Even more so, the value of is not 
measurable on temperature anisotropics alone. 

Similarly, tensor modes contribute a small amount of E polarisation. 
In principle, as the cosmic variance on E is smaller than on T, some 
constraints on T/S and/or can be obtained from the large scale 
measurement of Cf^. However, it is much safer to detect directly the 
presence of B modes of polarisation on large scales — an unambiguous 
imprint of inflationary gravity waves. As scalar perturbations do not 
produce odd parity polarisation {B modes), the detection of such polar- 
isation would be a smoking gun for primordial gravity waves, permitting 
to constrain the physics of inflation. 

There is, however, a caveat. It is quite possible that the ratio T/S is 
small enough to escape detectability (because of foreground polarisa- 
tion from lensing and foregrounds (see below), and of cosmic variance 
on the CMB temperature and E polarisation spectra. Still, polarisation 
measurement seem at present to be one's best bet for constraining 
inflationary models — even if the outcome is an upper limit on the 
energy scale of inflation. 

2.3.4. Lensing 

Lensing of the CMB polarisation field by large scale structure, which 
produces shear distorting the polarisation pattern, mixes the primor- 
dial E and B modes generated at last scattering. As a consequence, a 
contribution to B modes on small scales is expected even if there is no 
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primordial B. This effect permits to constrain the large scale structure 
power spectrum, quite interesting for constraining structure evolution 
through the "dark ages" at high redshift. 

Unfortunately, this effect also makes the detection of primordial B 
modes more difficult. If constraining inflationary B modes is the prime 
target of a CMB polarisation experiment, the cosmic shear imprint 
on the CMB polarisation is considered as a contaminant. Prior infor- 
mation on the large scale structure power spectrum, combined with a 
measurement of the E modes of polarisation, permits to predict and 
subtract statistically some of this lensing contribution to Cf^ (Knox 
& Song, 2002). 

3. Experimental issues 

3.1. Sensitivity 

One of the key issues for measuring CMB polarisation is the extreme 
sensitivity required, several orders of magnitude better than that re- 
quired to measure CMB anisotropics. 

At present, the best bolometric detectors built for sensitive balloon- 
borne CMB experiments are already photon noise limited. Therefore, 
reaching sensitivities better by two orders of magnitude (insufficient 
still for the ultimate polarisation experiment!) cannot be achieved by 
reducing detector noise. The gain in sensitivity requires about 10,000 
times more detectors or 10,000 times more integration time than present- 
day balloon-borne or ground-based experiments. There are essentially 
two strategies for doing so: either build a very large camera with 10,000 
detectors observing from the best possible site on ground, or multiply 
the number of detectors by an order of magnitude or two and benefit 
from the very long integration time that can be obtained on a space 
mission as compared to balloon experiments. 

Planck detectors will be somewhat more sensitive than present best 
balloon-borne detectors - by a factor of few - essentially because of 
reduced background and hence reduced photon noise. In addition, it 
will have about ten to twenty times as many useful detectors as present 
balloons - which essentially mesure anisotropics using the one or few 
best detectors on board. The total integration time will be about thirty 
times that of Boomerang and 1000 times that of Archeops. Although 
not all of Planck detectors are polarisation sensitive (the final mapping 
of temperature anisotropics is the primary goal of the mission), Planck 
will have the capacity to map the full sky polarisation of the CMB 
well enough for meaningful constraints to be put on both EE and BB 
polarisation spectra. 
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3.2. Foregrounds 

For CMB anisotropies, there exists a window in the electromagnetic 
spectrum for which CMB dominates over foregrounds on a large frac- 
tion of the sky. In addition, the near-independance of most emission 
laws on sky coordinates make component separation possible. 

The situation is far from being as clear for polarisation. At low 
frequencies, Faraday rotation complicates the model of the emission. 
At all frequencies, the integration of polarised emission of foregrounds 
along the line of sight results in some averaging of the polarisation, 
which depends on the spatial coherence of polarisation directions. In 
contrast, intensities add-up as long as the emitter is optically thin, 
which is mostly the case for CMB foregrounds. As a result, the polar- 
isation spatial power spectra of foregrounds can be different from the 
intensity spatial power spectra. 

In all cases, sensitive multifrequency observations are required for 
assessing the situation as far as foregrounds are concerned. Current 
experiments provide useful measurements of foreground polarisation 
properties. Recently, the Archeops team has observed the dust polar- 
isation at 350 GHz, finding 5-6% dust polarisation on average in the 
galactic plane, and as much as 20% in a few extended clouds (Benoit 
et al., 2003). 

Hence, because foreground polarisation monitoring is an important 
aspect of CMB polarisation measurement, it is important to measure 
polarisation at several frequencies. Such multi-frequency measurements 
allows one to check that the detected polarisation has the expected 
derivative of a blackbody emission law. If necessary, it also permits to 
evaluate and separate foreground polarised emission from the observa- 
tions. Planck will map the polarisation of the millimeter sky in six or 
seven frequency channels, which is not so easy to do from the ground 
(and depends on the polarisation properties of the atmosphere — not 
completely known yet). 

3.3. Multi-detector observations 

Presently, the most sensitive operational detectors for CMB polarisa- 
tion measurements are Polarisation-Sensitive Bolometers (PSB). These 
devices behave as a system of two polarimetric detectors, each of which 
is able to measure the intensity of the incoming radiation in one single 
linear direction. Each such polarimeter measures ideally a signal given 

by 

s{a) = 1-{I + Qcos2a + Usm2a) (10) 



delabrouille-polar.tex; 2/02/2008; 5:42; p. 9 



10 



Jacques Delabrouille 



where a is the angle between the polarimeter orientation (i.e. the 
polarisation sensitivity direction) and the axes chosen for the Stokes pa- 
rameters representation. For a PSB, the two polarimeters are oriented 
at 90° relative angle, and hence measure signals 

si = + Q cos 2^1 + [/sin2ai) (11) 

and 

S2 = -(/- Qcos2ai - C7sin2ai) (12) 

where a\ is the angle between the orientation of the first of the two 
polarimeters and the axes chosen for the Stokes parameters represen- 
tation. 

Whereas a single detector is sufficient - in principle -to make a 
map of CMB temperature anisotropics, it is not the case for measuring 
polarisation, unless the same polarimeter can measure polarisation in 
different directions at different times. A single PSB, consisting of two 
orthogonal polarimeters, cannot either measure both Q and [/ in a 
single observation. 

Whether one or several polarimeters are used in the end, measuring 
polarisation involves combining several measurements to separate the 
/, Q and IJ contributions. Clearly, at least three such measurements in 
each pixel are required to recover /, Q and JJ from the data. In general, 
a number m of independent measurements permit to recover the Stokes 
parameters by inverting the system: 

/ 51 \ 

Pure polarisation signals are obtained from differences of the read- 
outs of two polarimeters oriented at 90° , or possibly from more complex 
linear combinations of more polarimeter data streams. 

In a more compact form, recovering /, Q and ?7 at a given point p 
from all readouts (PSB and other) amounts to inverting a system of 
the type 

y(p) = ^5(p) + iV(p) (13) 

where Y is the vector of all measurements at point p, S is the vector of 
Stokes parameters at point p, 5 = (/, Q, [/), and N the vector of noise 
for the measurements. In the linear combination, lines of matrix A are 
ideally of the form Ai oc (1, cos 2q:^, sin lai) for a detector from a PSB, 
and Ai oc (1,0,0) for unpolarized detectors. 



/I 
1 



cos 20L\ 
COS 2^2 



sin 2^1 ^ 
sin 2^2 




\ 1 cos 2ara sin la^a / 



n2 
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The simplest approach to inverting the system is to use the pseudo- 
inverse = of A. If however, the estimated value A^ of 
A^ (computed from an assumed value A of A) is wrong, then A^A IS 
different of identity, and small off-diagonal terms induce a leakage of / 
into polarisation signals. The same is true (and worse) when one deals 
with T, and 5 as T > £; > 5 for CMB signals (at least in some 
of the useful regions of the spatial frequency range). Any inversion of 
the system with a matrix W (using the pseudo-inverse w = aK or 
the Generalised least square solution W = [A^N~^A]~^A^N~^^ or a 
Wiener inversion) suffers from the same problem of non-diagonal terms 
in WA when the estimate of 1^ is imperfect. This will be further 
discussed in the next section. 



4. The Planck design 

4.1. The Planck mission 

The Planck mission, to be launched by ESA in 2007, has been designed 
primarily for the ultimate mapping of CMB temperature anisotropics, 
with an angular resolution of about 4.5 arcminutes, and a sensitivity 
of about 2 X 10~^ per resolution element. 

The Planck optics comprise a 1.5 meter useful diameter off- axis 
gregorian telescope, at the focal plane of which are installed two com- 
plementary instruments. The HFI (High Frequency Instrument) is an 
array of bolometers cooled to 0.1 K with a dilution fridge, observing the 
sky in six frequency channels from 100 to 850 GHz, with polarisation 
sensitivity at 350, 220, 150, and possibly 100 GHz (decision pending 
for the last channel) . The LFI is an array of radiometers observing the 
sky at 30, 44 and 70 GHz with polarisation sensitivity in all channels. 

Planck will observe the sky from the L2 Sun-Earth Lagrange point, 
in a very stable thermal environment, away from sources of spurious 
radiation due to the Earth, the Moon and the Sun. The scanning is 
made along large circles around an anti-solar spin axis at a rate of 1 
rpm. The spin axis follows roughly the apparent motion of the Sun, so 
that the full sky is covered in slightly more than 6 months. 

4.2. Principles of polarisation measurement 

All Planck detectors, both in the High Frequency (bolometer) Instru- 
ment (HFI) and in the Low-Frequency (radiometer) Instrument (LFI) , 
are total-power detectors. 

In the LFI, all detectors are polarisation sensitive. In the HFI, some 
detectors are polarisation sensitive. The polarised detectors measure a 
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Figure 1. The two detectors in a single PSB share the same optics, but have different 
readouts. 

signal given by eq. 10, the other detectors measure a signal directly 
proportional to /. 

4.3. Instrumental set-up 

The Planck detector set-up takes advantage of the possibility to mea- 
sure two orthogonal polarisations using the same feed horn, by splitting 
the orthogonal polarisations into two different detectors. For each such 
feed, a pair of detectors sensitive to orthogonal polarisations (polarime- 
ters) share the same optics (filters, waveguides, corrugated horns, and 
telescope), but have different readouts as shown in figure 1. One single 
such device, hence, produces two different TOI, corresponding to mea- 
surements integrated in principle over the same beam shape, but with 
different polarisations. 

For the HFI, the splitting is done with Polarisation Sensitive Bolome- 
ters (PSB), in which two sensors located in the integration cavity are 
sensitive to orthogonal polarisations. 

As shown by Couchot et al. (1999), the measurement of polarisation 
using a set of polarimeters is most efficiently performed with so-called 
"Optimal Configurations" (OC), using a set of m polarimeters oriented 
at angles regularly spaced between and tt. A set of two feeds, with each 
two orthogonal polarimeters, oriented at 45° relative angle, provides 
such an OC. In the Planck focal plane, most (and possibly all in a final 
set-up) detectors are arranged in such a way that two complementary 
polarimeter pairs scan the sky on the same scan path, as illustrated in 
figure 2 
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Figure 2. The scanning of the polarised detectors permit to measure directly 
polarisation along four directions along each scan path. 

4.4. Aspects of data processing 

The reduction of Planck data involves several steps. The main data 
products to be obtained are: 

— calibrated maps of the sky (intensity and polarisation) in all the 
Planck frequency channels; 

— maps of the sky (intensity and polarisation) for cosmological and 
astrophysical diffuse emissions (CMB, thermal galactic dust emis- 
sion, synchrotron emission, ...); 

— catalogs of objects (clusters of galaxies, infrared galaxies, radio 
sources) ; 

— spatial power spectra of all components. 

Level 1 of Planck data reduction consists in putting all raw data 
(both scientific and housekeeping) into clean and organised Time Or- 
dered Information (TOI). Level 2 consists essentially in the map-making 
step which combines the level 1 outputs into calibrated maps for each 
frequency channels, as free as possible from noise and systematics. Level 
3 is the component-separation step, and level 4 consists in delivering 
the data to the community at large. 



delabrouille-polar.tex; 2/02/2008; 5:42; p. 13 



14 



Jacques Delabrouille 



We now discuss some of the polarisation-specific issues of Planck 
data-processing, with an emphasis on Planck HFI, although much of 
the discussion is appropriate for LFI data as well. 

Data processing involves three key aspects: modeling the problem, 
solving the problem, implementing the solution. Of these three aspects, 
we will now discuss to some extent the first two. 

4.4.1. Modeling the problem 

Modeling the problem consists, knowing the science issues and the 
way the data is measured (the instrument), to analyse the dependence 
function F in the relation: 

data = ^(instrument, science) (14) 

Both the instrument and the science can be modeled with a set of 
relevant parameters, some of which are known a priori, some of which 
are unknown. For some of the unknown parameters, statistical priors 
are available (e.g. uncertainties, allowed ranges, etc.). 

A single sample at time t from one detector in a PSB can be written 

as: 

y{t)=h{t)^s{t) + n{t) (15) 

where h(t) is the response of the detection chain (bolometer and read- 
out electronics), n{t) is noise including systematics, and s{t) the signal 
from the sky as observed through the optics: 

s{t) = J dndiy (ll{t) + QQ{t) + UU{t) + VV{t)) (16) 

In this equation, /, Q, [/, and V are the instrumental responses, as- 
sumed to be time-independent to first order, /(t), Q(t), U{t) and V{t) 
are the Stokes parameters^ o^th^ sky in the appropriate coordinate sys- 
tem. The dependence of /, Q, [/, V and /, Q, [/, V on sky coordinates 
is not written for readability. 

This model reaches a reasonable level of sophistication, although it 
is yet far from complete. Systematics are not specified; the response is 
assumed to be linear; sampling effect are not described in full detail 
(although the integrating window can be taken into account in the 
expression of h{t)). 

It is convenient to express all quantities in the instantaneous spin- 
axis frame, in which case a time delay 5t can be transformed in a 
coordinate shift as 50 = '0 = uJs\)\n^t- Then, the convolution h{t) s{t) 
becomes an integral on the sky, and we have for a stable spinning: 
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yit) = 



+n(i) 



(17) 

Where i ranges from to 3, So = Si = S2 = U and Ss = V and 
Stokes parameters are expressed in a sky reference frame for which the 
north pole is along the spin axis. 

The response to the sky signal for a stable spinning period (period 
for which the spin axis direction and the rotation vector magnitude 
do not change with time) is the polarized convolution, in the spin- 
axis frame, of the instrumental beam with the sky smoothed in the <p 
direction by the impulse response h. 

The connection of this equation with equation 10 is simple: Instead 
of being a known linear combination of Stokes parameters at a single 
point, measurements y are a weighted sum of Stokes parameters in 
a neighborhood in space and frequency (i.e. inverse wavelength). The 
responses Si and h{t) specify the elements of matrix A of equation 
13. Equation 13 is an approximation valid under the assumption that 
Stokes parameters are constant in that neighborhood. 

Polarinieters do not measure pure polarisation signals. The intensity 
response / is always larger than (or of the same order of magnitude as) 
Q and U for a polarimeter in a PSB, and the sky intensity is much larger 
than its polarisation. Neglecting the V term assumed to be negligible, 
combining polarimeter measurements into pure polarisation signals Q 
and/or U amounts to eliminating / from a set of equations as eq. 17, 
in which responses /, f/, V and h are not so well known a priori. 
This can be viewed as a component separation problem in which one 
of the components (the temperature signal /) is much stronger than 
the others and in which the mixing matrix is not perfectly well known. 

Hence, polarisation measurements simply "falling in the same pixel" 
cannot be straightforwardly combined to get polarisation by inverting 
the approximate model of equation 13, because of the risk of a signif- 
icant leakage of intensity / signals into polarisation signals (imperfect 
separation) in several ways: 



— If beams are different (absolute difference), the inversion yields an 
error proportional to / and to the beam difference. 

— If beams are identical, but oriented differently (case of a detector 
with an asymmetric beam "coming back" to the same point with a 
different orientation), the inversion yields an error proportional to 
/ and to the difference of the beams in the two orientations (which 
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has a quadrupole component, mimicking the effect of a quadrupole 
at last scattering). 

— If optical beams are identical and similarly oriented, but time re- 
sponses h{t) are different, the effective beams for the measurements 
are different. This yields an effect similar to the two previous cases. 

— If optical beams are identical and similarly oriented, time response 
effects assumed to be negligible, but pointings slightly different in- 
side the same pixel, then the inversion yields an error proportional 
to the gradient of / inside the pixel. 

— If the shape of the responses /, Q, [/, 1/ as a function of frequency is 
different (different spectral responses), the inversion yields an error 
proportional to / integrated in the spectral response difference. 

Note that although such effects exist also for temperature measure- 
ments, they are much more critical for polarisation, the effect being 
typically small as compared to the magnitude of /. Such effects that 
are below sample variance and foreground-induced errors for / may be 
above sample variance and foreground-induced errors for polarisation. 

4.4.2. Solving the problem 

Solving the problem amounts to finding the theoretical solution for 
estimating at best the useful parameters given the data and the priors. 

In Planck data processing, we take advantage from the fact that 
PSBs (or LFI radiometers) comprise two polarimeters sharing the same 
optics. Accurate models of the optical responses and measurements of 
the detector performance with Archeops and Boomerang indicate that 
the responses / of the two polarimeters in a same FSB are very similar, 
minimising the impact of the above systematics in the data as far as 
direct differences of the TOI from the two detectors in the same horn 
are concerned. 

Such differences are direct measurements of Q and U in the instan- 
taneous spinning reference frame (the frame for which the North pole 
is along the spin axis). These measurements can then be reprojected 
on polarisation maps directly, avoiding most of the / leakage problem. 

This nice design feature, however, does not solve all the problems of 
polarisation data processing. Some aspects of the data processing are 
discussed in Revenu et al. (2000), which investigate the suppression 
of stripes due to low frequency drifts, in Rosset et al. (2003), which 
investigate the effect of imperfect beams, of pointing inaccuracies, and 
of response mismatches in general, in Prunet et al. (2000), which dis- 
cuss a Wiener-based component separation of polarised emissions from 
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Figure 3. The four non-vanishing elements of the multivariate power spectrum of 
CMB temperature and polarisation for a particular cosmological model, with error 
boxes expected from Planck observations. 

CMB and foregrounds assuming prior knowledge of emission laws and 
power spectra, in Baccigalupi et al. (2003), which discuss the bhnd 
separation of polarised foregrounds specifically, in Lewis et al. (2002), 
which discuss the estimation of E and B spatial power spectra on 
incomplete sky maps. 

4.5. Science with Planck 

Predictions of the accuracy on temperature and polarisation spectra 
with Planck are shown in figure 3 (adapted from (Hu & Dodelson, 
2002)). Whereas error boxes are large especially on meaningful 
constraints can be put on cosmological scenarios. Note that this figure 
assumes an optimistic energy scale of inflation Ei — 2.2 x 10^^ GeV, but 
reionization later than what announced from recent WMAP results. 



5. Conclusion 

Measuring CMB polarisation is the next objective of the study of the 
cosmic microwave background. There are several motivations for doing 
so. 
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First, polarisation power spectra depend on cosmological parameters 
in a different way than the temperature power spectrum. Hence, they 
provide complementary information about cosmology to temperature 
anisotropics alone, lore and lore useful as the fundamental limitations 
come from cosmic variance and parameter degeneracies rather than 
experimental noise. In particular, early reionisation has a character- 
istic imprint on low I modes of CMB polarisation and temperature- 
polarisation power spectra. Also, polarisation permits to separate CMB 
anisotropy contributions induced by Doppler shifts (which do not pro- 
duce polarisation) from those intrinsic to density and tensor perturba- 
tions (which do produce polarisation). 

The B modes (odd parity) CMB polarisation field at low I can be 
generated solely by primordial gravity waves (in inflationary scenarios). 
At high ^, the weak shear due to structures in the Universe distorts the 
polarisation pattern, turning a small fraction of the E type polarisation 
into B type. Hence, measuring polarisation on all scales permits both 
to constrain the amount of gravity waves from inflation, as well as 
constrain the power spectrum P(k) of large scale structure, if the two 
effects can be properly disentangled (which is possible if T/S is not 
smaller than about 0.001). 

The Planck mission will map the polarised sky in six or seven fre- 
quency channels (30, 44, 70, 150, 220, 350 and possibly 100 GHz). The 
sensitivity of Planck to polarisation, of the order of ten microkelvin 
per resolution element on average and about ten times more in highly 
redundant patches, will permit to put strong constrains on the polari- 
sation power spectra of the CMB. 
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